Hyperfine interaction induced dephasing of coupled spin qubits in semiconductor 

double quantum dots 
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We investigate theoretically the hyperfine-induced dephasing of two-electron-spin states in a dou- 
ble quantum dot with a finite singlet-triplet splitting J. In particular, we derive an elTective pure 
dephasing Hamiltonian, which is valid when the hyperfine-induced mixing is suppressed due to the 
relatively large J and the external magnetic field. Using both a quantum theory based on resum- 
mation of ring diagrams and semiclassical methods, we identify the dominant dephasing processes 
in regimes defined by values of the external magnetic field, the singlet-triplet splitting, and inhomo- 
geneity in the total effective magnetic field. We address both free induction and Hahn echo decay 
of superposition of singlet and unpolarized triplet states (both cases are relevant for singlet-triplet 
qubits realized in double quantum dots). We also study hyperfine-induced exchange gate errors for 
two single-spin qubits. Results for III-V semiconductors as well as silicon-based quantum dots are 
presented. 

PACS numbers: 03.65.Yz, 76.30.-v, 71.70.Jp, 76.60.Lz 



I. INTRODUCTION 

Spin qubits in quantum dots (QDs) or donors have 
been extensively investigated during the past decadeiii"— 
Single spin qubits^ have been initialized, manipulated, 
and read out both electrically in gated dots^"^^ and op- 
tically in self-assembled dotsj^^"— Two-qubit exchange 
gates have been demonstrated recently in gated QDsr^iii 
so are two-electron^^ and two-hole^^ couplings in a pair 
of vertically stacked self-assembled QDs. Furthermore, 
a subspace formed by the singlet S and the unpolarized 
triplet To states of two electron spins has also been ex- 
plored as a logical qubit^^"— and coupling between two 
such qubits^ and their entanglement^^ has been demon- 
strated recently. 

Slow decoherence relative to the control speed is one of 
the key criteria for a scalable quantum information pro- 
cessor. In QDs made of III-V compounds or from natural 
Si, the dominant source of single-spin decoherence is the 
nuclear spin bath, coupled to the carrier spins by hyper- 
fine (hf) interaction^^"— The role of this hf coupling 
is clear: the energy splitting of an electron spin is af- 
fected by the fluctuating Overhauser field of the nuclear 
spins. The nuclear spin bath is quasi-static due to its 
dynamics being much slower than the dynamics of the 
electron spins. Therefore, the strongest effect of the hf 
interaction is an inhomogeneous broadening in the distri- 
bution of electron spin splittings due to the random but 
static orientation of the nuclear spins, causing a decay 
of the electron free induction signal on the timescale of 
T2 (ii^^ The slow nuclear spin dynamics within the bath 
causes homogeneous broadening, or pure dephasing of the 
electron spin qubit, which is measurable in a spin echo 
experiment. This nuclear spin dynamics is either due 
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or, at 



to pure hf interaction at lower applied fields 

higher fields, due to dipolar interaction between nuclear 

spins i^i^i^ 

For two uncoupled spins, spin product states such 
as Iti) and |4,t) are the two-spin eigenstates. If the 
two spins are initially prepared in a singlet state (S) 
(lti> - ||t>)/y2r^'^'^ the random Overhauser field 
would strongly mix the S and Tq states, and the timescale 
on which this mixing leads to decay of the measured sig- 
nal is T2*, the same as that for a single-spin qubiti^i^i^ 
Similar to the single spin case,* the application of a Hahn 
echo pulse sequence removes the influence of the quasi- 
static nuclear fluctuations, and reveals the much slower 
decoherence of two independent single spinsj^i^ 

For two coupled electrons in a uniform effective mag- 
netic field, the singlet and triplet states are the two- 
electron eigenstates. At a finite exchange splitting J in 
a double quantum dot (DQD), it has been shown^i^ 
that hf coupling is less effective in mixing the S and Tq 
states. The decoherence effect of nuclear spin pair-flips 
has also been analyzed4i Recently, decoherence in the 
S-Tq subspace has been investigated at finite J in exper- 
iments on GaAs^^ and InGaAs^"^ DQDs. In the former 
work the dominant role played by the charge noise was 
uncovered, while in the latter paper the effects of charge 
noise were minimized, and a lower bound on T2 time due 
to interaction with nuclei was obtained. 

In this article we systematically investigate the hf- 
induced dephasing of two-spin states at finite values of 
exchange splitting J. We focus on the case of dephas- 
ing within the S-Tq subspace, in which full control over 
the qubit state is possible because of the creation of a 
stationary Overhauser field gradient;^ or a gradient of 
magnetic field generated by a proximal nanomagneti^i^ 



We investigate both the hmit of uniform effective field 
and of the finite effective field gradient, and we find that 
with an increasing magnitude of the gradient there is 
a smooth transition from strongly suppressed decoher- 
ence to decoherence that is similar to the case of single 
spin qubits. We also study the hf-induced decay in an 
Hahn echo experiment, which is possible with a control- 
lable field gradient."*^ Lastly, we analyze the hf-induced 
exchange gate error, when the two spins in a DQD are 
treated as two single-spin qubits. 

Our theoretical approach is based on first obtain- 
ing an effective pure dephasing Hamiltonian via an 
appropriate canonical transformation of the full hf 
Hamiltoniaui^"— iiii^i^ The effective Hamiltonian H is 
diagonal in the basis of the relevant states, allowing de- 
phasing calculations for superpositions of these states. 
When dealing with terms in H that are of second order 
in the transverse Overhauser field, we use the ring dia- 
gram theory2i*2^i^ (RDT) . The terms that are of first or 
second order in longitudinal Overhauser field are treated 
classically. We also compare the RDT results with calcu- 
lations based on semiclassical averaging over the quasi- 
static transverse Overhauser fields, underlining its con- 
nection to the short-time RDT calculationsi^i^ 

It is important to point out that there are many other 
decoherence channels beyond hf coupling with nuclear 
spins for electrons in a DQD. For example, at finite J the 
orbital degree of freedom is not completely frozen, so that 
its fluctuations can lead to spin decoherence. In particu- 
lar, charge noise could be an important, even dominant, 
source of dephasing, as suggested by theory^i^"— and 
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Our results presented in this pa- 



per are most relevant in an experimental situation where 
charge noise is not dominant. However, despite the pos- 
sibly critical role of charge noise in DQDs, we would like 
to stress that this noise can in principle be removed to 
a large degree (e.g. by different designs of samples or 
the gate circuitry), while the presence of nuclear spins 
is unavoidable in III-V materials, since all isotopes of 
the III-V elements have finite nuclear spins (unlike Si, 
where isotopic purification could, in principle, suppress 
the nuclear- induced spin decoherence). Furthermore, it 
has recently been shown that S and To states of an opti- 
cally controlled self- assembled QD molecule can be ma- 
nipulated in a regime in which their splitting is insen- 
sitive (in the first order) to charge fluctuations)^ The 
existence of such "sweet spots" was also predicted^i^ 
in gated QDs. Besides charge noise, the different charge 
distributions for the singlet and triplet states also allow 
electron-phonon interaction to cause dephasing in the 
S-Tq subspace^ii^i^ although this dephasing channel 
should be weak in a double dot if J is not too large<^ 

In addition to pure dephasing, an S-Tq qubit can un- 
dergo longitudinal relaxation through phonon emission, 
which is allowed by spin-orbit coupling^"— and/or hy- 
perfine mixing between the statesi^ Such dissipative re- 
laxation processes lead to both transitions between S and 
To, and leakages out of the qubit subspace (to T± states. 



for example). The characteristic timescales are most of- 
ten much longer— than the timescales considered in this 
paper, although recent calculations suggest that in some 
parameter regimes the S'-Tq transitions in GaAs DQDs 
can occur at the microsecond timescale.^^ One can, how- 
ever, exploit the large anisotropy of the relaxation rates 
with respect to the direction of applied magnetic field, 
and suppress these processes by an appropriate choice of 
the in-plane B field directioui^ 

We mention that although the single qubit decoherence 
is often theoretically studied in the literature in various 
contexts, there are few concrete analyses of multi-qubit 
decoherence simply because the multiqubit decoherence 
problem is technically difficult due to the many possi- 
ble decoherence channels for the entangled system when 
even a few qubits are coupled together. Our current work 
demonstrates that the simplest entangled spin system, 
namely, a system of just two exchange-coupled electron 
spins weakly interacting with an environment of nuclear 
spins, is theoretically challenging, even though the full 
Hamiltonian for the problem is completely known. Real 
experimental situations are obviously far more complex 
because it is unlikely that all the environmental inffu- 
ences would be completely known. For example, as we 
mentioned above, for coupled spin qubits in semicon- 
ductor quantum dots, there would be, in addition to 
the nuclear-induced Overhauser noise, other decoherence 
mechanisms such as charge noise, stray fluctuating mag- 
netic flelds arising from random impurity spins in the 
semiconductor and from microwaves and other random 
fluctuations in the background. But the eventual con- 
struction of a fault-tolerant practical quantum computer 
necessarily requires understanding (and if possible, miti- 
gating) all multiqubit decoherence mechanisms since the 
quantum error correction threshold is small, which im- 
plies that only the smallest amount of decoherence can 
be efficiently eliminated by the error correction proto- 
cols. Our current work should be construed as a first 
step toward the goal of a comprehensive understanding 
of multiqubit decoherence in one of the most practical 
and widely studied proposed quantum computer archi- 
tectures, namely, the spin quantum computer in semi- 
conductor quantum dots. It is somewhat sobering that 
even this first step of understanding Overhauser noise in- 
duced two-qubit decoherence in semiconductor quantum 
dots is already a very challenging problem. 

The present paper is organized as follows. In Section 
im we describe the Hamiltonian of the two electrons in 
a DQD applicable in the regime of our interest. At fi- 
nite J, the Hamiltonian of two electrons is given in the 
{S', ro,r+,T_} basis in subsection III Al The state of the 
nuclear bath is described in subsection IIIBI In subsec- 
tion llTC] we derive an effective Hamiltonian in the S-Tq 
basis. The dephasing of an S-Tq qubit due to the various 
terms in the effective Hamiltonian is calculated for DQDs 
made of various materials in Section IIIIl At the end of 
this Section we identify the dominant dephasing mecha- 
nisms for various types of DQDs (GaAs, Si, or InGaAs) 



in subsection IIII El and discuss the functional character- 
istics of coherence decay that can be expected for various 
parameter regimes . In Section IIVI we clarify the effect 
of a finite interdot field gradient, and analyze the de- 
cay of the Hahn echo signal. The SWAP gate error due 
to the hf-induced dephasing processes is investigated in 
Section |Vl with a focus on the effect of inhomogeneous 
broadening. A description of experimental protocols to 
measure S-Tq decoherence is given in Appendix [^ Ad- 
ditional technical details are provided in Appendices B 
through D. 



II. TWO ELECTRON SPINS IN A DOUBLE 
QUANTUM DOT 

In this Section we define the starting point of our cal- 
culations. Specifically, we first identify the four relevant 
lowest-energy two-electron states in a DQD, then project 
the total Hamiltonian onto the basis that is a product of 
these four electronic states and the nuclear spin Hilbert 
space. We also discuss the semiclassical description of the 
nuclear spin reservoir. Starting from the general Hamil- 
tonian, depending on the particular physical problems, 
we derive the effective Hamiltonians that couple the two 
levels that we are interested in to the nuclear spin reser- 
voir. 



A. Low-Energy Tw^o-Electron Hamiltonian in a 
DQD 

The system we consider in the current study consists 
of two electrons located in two (weakly) tunnel-coupled 
QDs, labeled left (L) and right (R), deep in the (1,1) 
charge configuration . ^^■'^^ —'^^ The total Hamiltonian for 
the coupled electron-nuclear spin system can be written 



H, 



total 



^H, 



H, 



hf 



Hn ; 



(1) 



where the three terms represent the electronic, hyperfine, 
and nuclear part of the Hamiltonian, respectively. Below 
we discuss each of these terms in a two-electron DQD. 

Deep in the (1,1) regime of the charge stability dia- 
gram, the four lowest-energy two-electron states, includ- 
ing one singlet 5(1, 1) and three triplet states T{1, 1), are 
well approximated by the Heitler-London states. Specif- 
ically, 5(1, 1) = Vs ® (lU) - lit))/v^ and T(l, 1) = 
V'AS «) Itt) , (lU) + \it))/V^, and lU)- The orbital parts 
V'S/AS are symmetric and antisymmetric combinations of 
^L(r) and ^_R(r) states, which are the single-electron 
ground state orbitals of the potentials for the L and R 
dots. States in the doubly occupied (2, 0) and (0, 2) con- 
figurations have significantly higher energies, and are not 
included explicitly in our consideration. The only role 
they play is to lower the energy of the singlet with re- 
spect to To by the exchange splitting J (J>0 for typical 



values of magnetic fields used in experiments) , when the 
DQD has a finite tunnel coupling j^^'^^ 

In the presence of a magnetic field, which has both 
a uniform B — {Bl + Bb)/'2, and a gradient Ai? = 
[Bl — B^)l2 component, the electronic Hamiltonian in 
the {S',To,T+,r_} basis is 



-ffp = 



-J b.Ez 

l\Ez 

~Ez 

£;z 



(2) 



where Ez — —gfisB and AEz = —guB^B (the sign 
convention is such that the positive B field lowers the 
energy of T+ in GaAs, where the effective g- factor is neg- 
ative). 

The two electrons couple to the environmental nuclear 
spins through the contact hf interaction, which takes the 
general form 

Hhf =^Aaii]Si-Iii^oS{ri-'Ri)+^Aa\i]S2-IiiyoS{r2-Iii) , 

i i 

(3) 
where Si, 2 are the spin operators of the two electrons 
at positions ri^2, h are the spin operators of nuclei at 
site Hi, and Aa[i] is the hf constant corresponding to the 
species a[i] (e.g. ^^Ga, '^^Ga, and '^^As in GaAs, or ^^Si 
in Si) of the nucleus at site i. The values of hf constants 
for relevant nuclei are given in Table HI i/g is the vol- 
ume of the primitive unit cell, and using a single-electron 
wavefunction 0(r) = — ^vl/(r), where ^(r) is the enve- 
lope function, the hf interaction energy (i.e. the Knight 
shift) of the i-th nucleus interacting with one electron is 
A, =^„[,j|^(r,)|2. 

Projecting Hamiltonian ([3]) onto the {S,To,T+,T-} 
basis (see Appendix [B] for details), we obtai n'^^i^°i^^ 



Hhf - 
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(4) 

(5) 
(6) 



'I'i(j^)(Rj)P denotes the hf coupling 



Here A\^^^ = A^ 
between an electron in the E{¥C) orbital and a nuclear 
spin at Ri . Here we have neglected the overlap between 
the L and R orbital wavefunctions. Keeping the finite 
overlaps amounts to small quantitative corrections to the 
matrix elements of i?^,/, and neglecting them does not 
cause any qualitative change (for the general form of B^f 
and its derivation, see Appendix [Bjl . 

The diagonal terms for the two polarized triplet states 
in Bhf are the longitudinal Overhauser field along the 



external field direction z. On the time scale of interest 
to the present study, it is quasi-static. ?M0 We therefore 
express the field operator in terms of its ensemble average 
and fluctuations: 

A = ^ Cdl = Mo + -5/* (7) 

i 
i 

where (...) denotes an average over the nuclear bath. 

In Hamiltonian -ff^,/, the Overhauser field difference 
between the two dots, ^^ BJl, couples S and To states. 
It can be an important control for universal manipulation 
of an S'-To qubit;^ Again we split this term into the 
average and the fluctuations^ 



E^^^ 



(9) 



Here the mean field 6*0 = (6*) can be built up through 
dynamical nuclear spin polarization (DNP),^^ while the 
fluctuation S9 = 6 — 60 can be reduced with respext to 
its "natural" high-temperature value during the DNP 
process^ 

In experiments where the Overhauser field in the DQD 
are prepared by DNP through multiple sweeps across the 
S-T+ anticrossing)22i2^iM both finite 6*0 and mo are es- 
tablished. Typically both 9o and mo are of the order 
100 mT (or 2.5 MeV for GaAs). In Ref. IH 0o reaches 
above 100 mT, with no reaching approximately 130 mT. 
Note that an equivalent role can be played by an exter- 
nal magnetic field gradient AB^ (leading to finite AEz), 
which could be established using a nanomagnet located 
close to the DQD,^iiii^ with reported values of - 10 mT 
field difference between the two dots. 

The total Hamiltonian for the two electron spins and 
the hyperfine interaction now takes the form 
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(10) 
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Here we have combined the quasi-static mean-field Over- 
hauser terms with the external magnetic field, with AEz 
replaced by 



and Ez replaced by 



Mt = 



.E,_ 



AE, 



Mo 



(11) 



(12) 



TABLE I: The hf constants and the nuclear Zeeman energies 
for i? = 1 T. The parameters for Ga, As and In are taken 
from Ref. |33 . For ^^Si, the value of Ac comes from Ref. iGq . 
For comparison, —Ez/g ~ 57.8 ^eV at B = 1 T. 



Nuclear species a 



Ac (msV) 



(neV) 



«3Ga 
^iQa 
^5 As 
ii^In 

29 r 



^'Si 



35.9 


-42.1 


45.9 


-53.6 


42.9 


-30.1 


55.8 


-38.4 


56.0 


-38.5 


2.15 


34.8 



The last term in Htotai is the nuclear Zeeman energy: 



Hn — 2_^^a[i\Ii , 



(13) 



where uja^i^ is the Zeeman splitting of the nucleus of 
species a at site i. Typically these splittings are smaller 
than the electronic one by 3 orders of magnitude (see Ta- 
ble m . Note that including the finite values of Wa in the 
case of multiple isotopes (as is the case in III-V based 
QDs) is crucial for description of Hahn echo decay of a 
single electron spin^"— (or an S'-To qubit"^^ at J — 0, 
which is equivalent to two independent single spins), 
while the nuclear Zeeman energies generally have much 
smaller influence on dephasing during the free evolution 
of the qubit. Below we will show that these statements 
also hold for the S'-To decoherence at finite J. 



B. The nuclear bath and its semiclassical 
description 

As we have discussed at the end of the previous subsec- 
tion, one of the key features of a coupled electron- nuclear- 
spin system is the smallness of the intrinsic nuclear en- 
ergy scales (both the Zeeman energies and the dipolar 
interactions among the nuclei). Consequently, at experi- 
mentally realistic temperatures, the equilibrium nuclear 
density operator is proportional to unity, /5/ ex 1. When 
nuclear spins are dynamically polarized^^i^i^i^i^"— 
the direction of the polarization in each dot is along the 
applied field (z) direction. The components of the nu- 
clear spins transverse to this direction are randomized 
on a timescale of ^ 100 ms due to intra- nuclear dipolar 
interactions,--- so that for experiments in which the total 
data acquisition time is much longer than this timescale, 
the appropriate nuclear density matrix is diagonal in the 
basis of eigenstates of If. A semiclassical description 
of the nuclear reservoir is thus valid for at least some 
situations.— Here we discuss some of the most impor- 
tant characteristics of the nuclear reservoir. 

The Overhauser field is defined as h = '^^Aili. The 
maximal value of the Overhauser field (on a single elec- 
tron in a given orbital) in a fully polarized nuclear bath 



Am = ^ Aili = 22 "n-alaAa , 



(14) 



where i denotes the nuclear sites, It is the i-th nu- 
clear spin, a denotes the nuclear species, and n^ is 
the average number of nuclei of this species in the unit 
cell (i.e. in both III-V compounds and in Si we have 
X)q ""c ~ ■^)' ^^'^ -^a Eire the hf couplings of nuclei of 
a species given in Table HI With the envelope functions 
^L,B normalized as / |^L__R(r)pd'^r = j/q, with vq being 
the volume of the Wigner-Seitz unit cell, we have then 
A ' =-^a[j]l'^i,fl(i"j)Pj 9,s stated before. We also define 
the number of unit cells A^* in which the probability of 
finding an electron described by wavefunction ^(r) has 
appreciable magnitude 



iV* 



This definition implies that 
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(15) 



(16) 



where the sum over u is over all the Wigner-Seitz unit 
cells (where we assume that the envelope function is prac- 
tically constant within each cell). 

On the timescale when the nuclear spins can be con- 
sidered static, we can replace the quantum averages over 
the nuclear bath by classical averages over the values of 
the static Overhauser field h described by a Gaussian 
probability distribution)^ 



Tr 



PI f{Y,Ad^)]^ J P{h)f ih)<f 



(17) 



where 



F(h) = 



27rcrj_ 



^hl/2al_^_ 



-{h.~{h,)fl1al^ 



2'ITCTz 



(18) 

where hj_ is the transverse component of h, and {h^) is 
the average value of the longitudinal Overhauser field (in 
the z direction). The width of the distribution of h^ is 
given by 



'z,* 



^y,^.wt?)-{I^f 



(19) 



In the following we will simply replace {{If)'^) — (Zf )^ by 
its zero-polarization value [which, by isotropy, is equal to 
Ia{Ia + l)/3 for i G a] multiplied by a narrowin g^^'^^" — 
factor njp < 1 , which takes into account the narrowing of 
the distribution caused both by creation of finite polar- 
ization, and by processes which diminish the spread of 
the values of the Overhauser field without affecting its 
average value. Thus 
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where we have defined the "typical" energy of hf interac- 
tion A. For GaAs A is of the same order of magnitude 
as the maximal Overhauser field Am- For silicon with a 
fraction / of the spinful ^^Si nuclei, AocAm/VJ- 

For typical achievable values of nuclear polarization, 
the width of the distribution of the transverse Overhauser 
field is given by an analogous formula, albeit without the 
np factor. In other words. 



c^,* 



np 



(21) 



In III-V compounds, the number of nuclei interacting 
appreciably with an electron is Ns = 2A''^ (the factor 
of 2 appears because there are two nuclei per unit cell) , 
and the typical value a± for a GaAs QD with Ns ~ 10^ 
spins is a few mT (< O.lfieV). For silicon, the result 
depends also on the concentration of the spin-1/2 ^^Si 
nuclei, given by / = nsi/2 (with / = 0.047 for natural 
silicon). With /si = l/2, we obtain 



,Si 



*(/) - V</ 



-4s 



Asif 



^/2N^ 



(22) 



where the value of Asi is given in Table U (note that we 
are using here a different definition for ^si compared to 
Ref. 66, where an /-dependent quantity was used). 

The difference of the Overhauser fields between the 
two dots is now: 



(23) 



which could have a finite average and/or a narrowed dis- 
tribution. Using the values of <Jz,l/r ^ot the two dots (L 
and R), we introduce 

A^ 



where 



2 /'„2 

fK^z.l 
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in)- 
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Nd 
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1 



N 



(25) 



and we have taken the natural (non-narrowed) values for 
'^z.L/R- The average transverse fiuctuation a± can be 
similarly defined as Eq. ([24|) . but without the np factor. 



C. Effective Hamiltonian in the S-To subspace 

One focus of the present paper is the decoherence of 
5- To qubits. Here we derive the effective Hamiltonian 
in the basis of \S) and |To) states in the presence of a 
large external magnetic field. It is directly applicable to 
experiments on S-Tq qubits whenever exchange splitting 
J is large enough. 

As shown in Eq. ([T0|) . |5) and jTo) states are coupled to 
the polarized triplet states \T±) by the transverse Over- 
hauser field X^i ^i ^i • In a finite external magnetic 



field, such that /ix ^ crj_, dephasing between S and Tq 
states can be faithfully described by an effective Hamilto- 
nian in the subspace of these two states, treating the cou- 
pling to the T± states perturbatively. With the zeroth- 
order Hamiltonian given by the first matrix in Eq. (jlOp. 
the condition for the perturbative treatment is 



cr_L < |-^±Mt|, ImtI 



(26) 



The effective Hamiltonian in the {S,To} subspace is 
then>ii 



H 



STo 



-J + Vss Vst„+Ot 



K 



STo + ''T 



-"^^alill! 








(27) 



which contains second-order effective interactions among 
the nuclei. In particular, Vss comes from the virtual 
flip-flops between S and T±: 



^^^ = j^-TT^ 



J2BrB,I+I- 



= V,, J2i^^^^ + ^f^f - AlAf - AfAf)I+ir , (28) 

with Vgg = J/A{fM^ ~ J"^)- ^STo represents flip-flops be- 
tween S and To via virtual transitions through T±. It 
consists of a hermitian Vh and an anti-hermitian Vah 
part: 



VsTo —Vh + Vah , 



(29) 



where 



Mt 



4 \/iT fl^ ~ J^ 



J2b.c,{i+i- + i-i+) 
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VH^Zi'^^^f - AfAf)iltlJ + 17 It) (30) 

1 J 



^^" - "4m^-J2 



Y^B^CAltlJ-I^lt 



«j 



VAHY.^A^^f-AfA'^)ltr, 



(31) 



with 



1 / 1 



Mt 



16 \^J.T Mt ~ "i^ 



1 J 

VAH = ---2— 



Ai|- J2 



Figure [T] gives a cartoon that depicts the virtual flip-flop 
transitions contributing to Vss and Vstq ■ 

Hamiltonian Hsto in Eq. ([ST)) contains terms that are 
linear in the nuclear spin operators /f but are second- 
order in the hyperfine coupling strength (within our ap- 
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FIG. 1: The virtual transitions that lead to the second-order 
effective interactions Vss and Vstq, as described in text. 



proximation of neglecting the orbital overlap when cal- 
culating the hf interactions, Bf = Cf = A^/A): 
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4(AiT + J) 



7 , A^ /j , 



LE^^ = i^E^^^- 



Mt 



(32) 



They influence the S-Tq coherence in a way identical 
to what the longitudinal Overhauser fleld does to a sin- 
gle spin. However, the normal Overhauser field ^^ Ail^ 
leads to a strong inhomogeneous broadening^^ and a de- 
phasing time of Tj* ~ l/cz, where Uz is the spread of 
the values of the longitudinal Overhauser fleld (see Sec- 
tion [HB] for precise definition). This time is of the order 
of 10 ns in GaAs QDs for a single electron spin, while 
dephasing due to these /^-linear terms is strongly sup- 
pressed in the case of an S-T^ qubit, because the interac- 
tion strength is significantly reduced (i.e. the qubit-nuclei 
couplings in Eq. ([5^ are ~ Af/^T^Ai for the values of 
/iT considered here). 

Note that Vss, Vstq, and 69 are all fluctuations due to 
nuclear spins, and average to zero in a thermal nuclear 
spin reservoir. In the absence of these fluctuations, the 
universal control of the S-Tq qubit can be achieved via 
tuning of J and 6t-^^ In the context of decoherence when 
J is flnite (more precisely, when y^,P + 46'| ^ crz,a±), 
we have two interesting regimes to consider. The first is 
when 6't = 0, i.e. the DQD is in a uniform total field. In 
this limit, the DQD has left-right symmetry, so that the 
singlet and triplet states are system eigenstates. Since S 
and To are Sz = states, with no magnetic moment in 
either quantum dots, they are not directly affected by the 
Overhauser field hz ■ The S'-To qubit made up from these 
two states should thus have a significantly longer inho- 
mogeneous broadening time (T| ) than a single spin. We 
will study this regime in detail in Section Hill The other 
regime is when 9^ 3> <7z, i.e. a magnetic field gradient 
is present (whether due to nuclear spin polarization or 
applied externally) . Now the left-right symmetry of the 
DQD is broken, and the true eigenstates of the system 
are superpositions of S and To states. The electron spin 
densities in the two dots do not vanish anymore, so that 
hz can affect the two-spin coherence directly, and the 



system acquires single-spin qubit characteristics. This 
regime wih be studied in Section IIVI 

An extreme case is when J <^az, <J± , which is a regime 
aheady investigated in existing experimental2ii2^i2^i^i2S 
and theoretical^iil^ studies. Here the exchange couphng 
is effectively turned off. The dynamics of the two inde- 
pendent electron spins are determined by the Overhauser 
fields in the respective dots. The inter-dot nuclear spin 
fiip-flops are completely suppressed, so that all interac- 
tions involving such fiip-flops, Vss and Vah, vanish in 
this limit. In addition, now Hb ~ He so that they do not 
affect the two-spin dynamics. The remaining hf-mediated 
interaction is due to intra-dot fiip-flops: 



five magnetic field. 



lim Vh 

J-i-O 



'^ E(^'^' - A?Af){I+lr + I-ip . 



«j 



(33) 

When J = 0, the two-spin eigenstates are spin prod- 
uct states. Within the Sz = Q subspace, it is more 
convenient to consider the Hamiltonian in the basis of 
|±^) = ^(|5) ± |To)) = {| ti>, I it)} states. The result- 
ing Hamiltonian is of pure dephasing form in this product 
basis: 



i?j=o « [VH + eT+Sm+X) {+X\-\^X) {-X\) . (34) 

In a free evolution experiment with an initial singlet 
state, the ensemble coherence will now decay in T2* « l/cr^ 
due to the 59 term i^^i^^ On the other hand, in a Hahn 
echo experiment )Sii^ the influence of 59 is removed, and 
the signal decay is due to Vr from Eq. (1551) . Since this 
interaction is a sum of two commuting terms from two un- 
coupled dots, the appropriately defined S-Tq decoherence 
function is a product of the two single-dot decoherence 
functions.™^ This observation establishes a one-to-one 
correspondence between single-spin Hahn echo decay^ 
due to hf-mediated interactions considered theoretically 
in Refs. l30l - [33l47l and the J = singlet-triplet Hahn echo 
decayi^i^ 



III. CALCULATIONS OF S-To DEPHASING IN 

THE ABSENCE OF INTERDOT MAGNETIC 

FIELD GRADIENT 



In this Section we study two-spin decoherence within 
the S-Tq subspace in the absence of interdot magnetic 
field gradient, i.e. 9^ =0 (practically, this is true as long 
as 0T < CTz)- With J » cr^, cr_L (note that the magnitudes 
of Vss and Vstq are negligible compared to Uz , as shown 
in Appendix[C]) , we can perform a second canonical trans- 
formation to diagonalize Hamiltonian (j27p . treating the 
off-diagonal terms as a perturbation. We now obtain the 
final form of the effective Hamiltonian in a uniform effec- 



Hg,j.^Q — 



-J + Vss + HA + Vse 







-Ha - Vse 



-2_/^"H'^^ 




where 



Ha - ~^T.B^B,m] = -^ 

i-'j 

Vse = -J {{Vh, 59} + [Vah, 59] 



(35) 



(36) 



(37) 



To evaluate Vgg, we use a 1/N type approximation where 
we neglect the commutators of nuclear operators. The 
commutator in the above equation thus vanishes, so that 



Vse ~ ve I 



where 



ij 



Ve 



{AfA 



AfAf){I+r 



fJ-T 
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(38) 



8 J Vmt /^t ^ J^ 



With 0T = 0, S* and Tq are two-spin eigenstates. The 
dynamics of the S'-Tq coherence is quantified by the fol- 
lowing decoherence function 



WstM 



_ I aiQ 



it) 



STa 



(0) 



= Tr, 



{pie'"' 



Tot^-iHst 



(39) 



Here /5/(0) is the initial nuclear density operator, while 
Hs and Htq are the operators appearing on the diagonal 
in the effective Hamiltonian from Eq. (j35l) . In Appendix 
El we discuss how the above quantity can be measured in 
electrically controlled S-Tq qubits. Some of the informa- 
tion contained in Wstq [t) function may also be indirectly 
inferred from the optical spectra of two-electron states in 
coupled self-assembled InGaAs quantum dots4^ 

Our main goal is to quantify each of the dephasing 
processes as we vary J and /it- There are three kinds 
of hf-related terms appearing in Herj.=Q'. (1) Hb and He, 
the /^-linear terms, (2) Ha, the square of the longitudinal 
Overhauser fleld, and (3) Vss and Vse, which are second 
order in transverse Overhauser field. Below we discuss 
their individual contributions to S-Tq dephasing dynam- 
ics. Such a treatment makes sense when the timescales 
on which various terms operate are very different. If the 
importance of these interactions is comparable, treating 
various interactions as independent would introduce a 
quantitative error. In the following we will present re- 
sults for the envelope W{t)^ in which we have removed 
the fast oscillating part of WsTa (*) from Eq. 



Wit) 



= e-^^* 



WsToit) 



(40) 



A. Dephasing due to Hb and He 

While Hb ~ He leads to dephasing that is completely 
analogous to the hf-induced inhomogeneous broadening 
of a single spin free evolution, the magnitude of the de- 
phasing here is strongly suppressed because of the re- 
duced coupling between the nuclei and the qubit. From 
Eq. (p9|) . keeping only Hb and He in the Hamiltonian, 
we obtain 



W^B.cW=Tri(p,(0)e'(^^-^«)* 



Neglecting the wavefunction overlap, He 

7(E.eL^f^f + E,ei?^?^f)-with 
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J 



4/zt(/^t + J) 



(41) 



Hn = 



(42) 



For a large number of nuclei {Ns ^ 1), as discussed in 
Section III B[ tracing over the nuclear spin density matrix 
in Eq. (|41|) can be approximated by averaging over a 
classical Gaussian distribution of energy splittings with 
variance 



'B,C 
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where the variance for the i-th spin, af — {{liY) — {li)'^, 
is approximated by its value of ■^lil + 1) at vanishing 
nuclear polarization p (this is a good approximation at 
small p, since corrections are of the order p^). The value 
of ^^ Aj depends on the distribution of the nuclear cou- 
plings, i.e. the shape of electron wavefunction, but it can 
be roughly estimated as 1/iV^, as shown in Appendix IdI 
The resulting coherence decay is then given by 



with 



Tbc 



WBc{t) = e 



V6 



_ ^-{t/T;? 



(44) 



\l\nFyJY^' 



^a^a \-^a 



l)Ai 



Nl' + N^' 



(45) 

where Nl and Nji are the numbers of unit cells in the L 
and R dots (as defined in Eq. ([T5|)). The typical values 
of this time for GaAs and Si are shown in Fig. [2] 

The very long Tbc shown in Fig. [2] is a clear illus- 
tration of the strongly suppressed electron-nuclear spin 
interaction, which is a consequence of the highly symmet- 
ric nature of the singlet and unpolarized triplet states. In 
the following we generally neglect the contributions from 
Hb and He- 



B. Dephasing due to the second-order longitudinal 
Overhauser field Ha 



ffyi-induced dephasing can be calculated classically be- 
cause of the quasi-static nature of the longitudinal Over- 
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FIG. 2: The timescale Tbc of S-To coherence decay due to 
J^-linear terms (see Eq. (|^) for (a) GaAs DQD with iVi, = 
Affl = 10*^ at B = 200 mT (|^t| ~ 5.0 ^eV), and (b) DQD 
based on natural silicon with Nl = A''^; = 10^ at B = 30 mT 
{\^J.T\ ~ 3.5 /xeV). The regions in which a± <^ \J i: /^tI, I^t 
is not fulfilled are hatched. 



hauser field. We can rewrite Ha in terms of the clas- 
sical Overhauser field 9 = {hf^ — ft,|j.)/2, where hi^/p = 



E.^ 



L/R-, 



Ha^--. 



(46) 



As in Section [II Bl we treat as a Gaussian random vari- 
able, and obtain the relevant decoherence function Wa {t) 
by evaluating a Gaussian integral: 



where 



WA{t)= p{e)e^''^'''Ue 



P{0) 



2Tr(Tg 



(47) 



(48) 



with a distribution width ae — cJz/'^, where CTz is defined 
in Eq. ([M]). We then obtain (see also Ref. lill) 



WAit) = 



62 



£irctan{rjAt) 



a+vit^) 



1/4 



(49) 



where we have defined rjA — 4:<^g/J — <^1/J- The 
characteristic decay timescale Ta{Ot^ = 0) is defined by 
\Wa{Ta)\ = 1/e: 



Ta = 



^2 






(50) 



In Fig.|3]we show the typical values of this dephasing time 
for GaAs and Si. These Ta dephasing times are much 
shorter than Tbc from the previous subsection, but are 
significantly longer than the single-spin inhomogeneous 
broadening dephasing time Tj* , which are in the order of 
10 ns in a GaAs QD. As such, Ta represents the inhomo- 
geneous broadening of the pure S-Tq two-level system at 



finite J when 9^ = 0- Qualitatively, Ta is long because 
Hji is second order in the Overhauser field fluctuations 
S6, while for single spins the inhomogeneous broadening 
is in the first order of S9. The very long Ta, as compared 
with the short single-spin Tg*, means that narrowing of 
the nuclear bath may not be necessary for S-Tq qubits if 
9t = 0, though it is worth noting that Ta does increase 
as rip when the nuclear bath is narrowed. 
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FIG. 3: The timescale Ta of S-Tq coherence decay, (a) GaAs 
DQD with Nl=Nr^ 10^ at 5 = 200 mT (|/iT| ~ 5.0 ^eV), 
and (b) natural Si DQD with Nl=Nr^ 10^ at B = 30 mT 
{\ht\ ~ 3.5 fieV). In both cases the effective field gradient 
8t~0- The hatched regions are where crj_ ^ [ J± /iT|, Ia^tI is 
not fulfilled. 



where Uk are nuclear Zeeman energy of the fcth nu- 
cleus. The Ha, Hb, and He terms would have modi- 
fied this time dependence by changing the nuclear spin 
frequency (in particular, keeping the nonlinear Ha term 
would make the exact calculations much more compli- 
cated). However, it turns out that the timescale on which 
the dependence of I^ (t) becomes visibly modified by the 
presence of one of these terms is typically much longer 
than the timescale of coherence decay due to either this 
term by itself or due to Vss- For example, including 
Hb and He in the interaction picture would lead to cor- 
rections in time-dependence of transverse nuclear oper- 
ators on timescale of is ~4iV|)(^T + J)/ J?, which will 
turn out to be much larger than Tss (derived below) 
unless the number of nuclei is very small (which could 
happen for example, in small QDs made of isotopically 
purified silicon) . Similarly, one can quickly estimate that 
the timescale on which the corrections due to Ha appear 
is much larger than Ta given in Eqs. (|50p . unless Ns or 
np are so small that np^/Ns < 1. In short, neglecting 
the /^-dependent hf terms does not lead to any qualita- 
tive change in our calculations below, while it dramati- 
cally reduces the complexity (and therefore increases the 
transparency) of the calculation. 

We now outline our calculation of dephasing due to 
Vss- The decoherence function is 



Wss{t) = Trj {piT 



-if^dTVss(r) 



(52) 



C. Dephasing due to Vss 

Ha, Hb, and He lead to inhomogeneous broadening 
because longitudinal Overhauser fields are quasi-static. 
On the other hand, Vss smd Vse are due to transverse 
Overhauser fields, and their influence cannot be fully re- 
moved by a Hahn echo (see Section llVBI) . In this subsec- 
tion we study S-Tq dephasing induced by V^s- 

According to our general approach, in which we treat 
each of the dephasing-inducing terms in the full Hamil- 
tonian separately, we write the effective Hamiltonian as 
H = (Hz + Vss) \S) {S\+Hz \To) {To\. The Zeeman term 
is kept here, because while it does not lead to any dephas- 
ing by itself, it modifles the dynamics caused by Vss, as 
we will show below. 

An important question is why we can neglect Ha , Hb , 
and He terms in the current calculation. These terms, 
apart from leading to S-Tq dephasing in ways described 
above, in principle also modify dephasing caused by Vss 
by giving different Knight shifts to nuclear spins at dif- 
ferent locations in a quantum dot. By neglecting these 
terms, the operators in Vss are expressed in the inter- 
action picture only with respect to the "noninteracting" 
part of the Hamiltonian. Speciflcally, the nuclear opera- 
tors in the interaction picture are: 



Ik it) 






(51) 



where T is the time-ordering operator, and Vss{t) is Vss 
in the interaction picture with respect to Hz, so that 
the nuclear spin operators are given by Eq. (|51l) . The 
calculation is performed using the l/A^-expansion-based 
RDT from Refs. 31,32. In the interaction picture. 



Vssit) = Vss Y^ (t)kAAkAiJ+{t)J-{t) 



(53) 



kd 



where 0i = 1 for z e L and (/),i = — 1 for i e i?. When aver- 
aging over the nuclear operators, we contract the nuclear 
indices in pairs, so that the sign factors drop out {(fyf = 1). 
The T-matrix defined in Ref. 32, after taking into account 
the slightly different structure of the averaged exponent, 
is then given by 



Ti 



kl 



s2y^akaiAkAi e 



iUkit/2 



t^kl t 



^kl 



(54) 



where Uki ^ujk-uji, and Uk = (Ikh >o = ihih + 1) for 



an unpolarized nuclear spin bath. Applying the RDT, 
we re-sum the linked-ring terms from the perturbative 
expansion and obtain a general formula^ 



wssit) = n 



i^ p — i ai'ctan(ATn (t)) 



y/i + {Xm{t)y 



(55) 



where Xm{t) are the eigenvalues of Tki. As discussed in 
Ref. |32| . we can simplify the calculation by introducing 
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a "coarse-grained" T-matrix of Nj x Nj dimension, with 
Nj being the number of distinct nuclear species a, each 
of them having a distinct value of the Zeeman splitting 



_i/lo/vA/; 



sm ■ 



ra/3 = 2w,,VJw^Vn^n^^^ e-=/it/2 

(56) 
For a system with multiple species in the short-time 
limit {t<^ l/iOafs, when Tap ^Sap), or for a system with 
a single nuclear species, we obtain 



Wss{t<t:i/LOap) 



p — i arctan(?7^^ t) 



(57) 



where 



IvssliY^akAl + y^akAl 



7 , ^K^'-k ' /_^ 

keL k£R 

2\vss\{cr]_j^ + alj^) = 2\vss\<j'i 



(58) 



The result of Eq. ([57)) can also be obtained using a semi- 
classical quasi-static bath approximation, as shown in 
Appendix [Cl If the characteristic decay time T^^ defined 
by |l^ss(rgs)| = 1/e falls in this short-time regime, we 
obtain 



Tss — 



Vi2 



2v, 






For a system with multiple nuclear species in the long- 
time limit, t 3> 1/wq^, the heteronuclear spin contri- 
butions in Eq. ([54)) become negligible compared to the 
homo- nuclear ones. Then Wss{t) can be approximated 
by a product of functions describing decoherence caused 
by each spin species treated separately: 



Wss{t:>i/ujap) ~ Y[WssAt) 



(60) 



where the homonuclear contribution with the nuclear 
species a is given by 



W^SS,a(i>l/Wa/3) 



g-iarctan())ss,Q*) 
^1 + iVSS.at)^ 



(61) 



with ?7ss.c( = 2\vss\naA.'^/ND. In GaAs, where all rja 
have similar values, the estimate of the characteristic de- 
cay time (when it indeed falls in the t^ l/ujap regime) 
is Tss ~ l/f^a- The asymptotic decay of the coherence 
function (for i^maxajy"^^) depends then on the num- 
ber of nuclear species, e.g. Wgg^{t) ^ t^'^ in GaAs, while 

WsT(t) ^ i"' in Si. 

In Fig. ID we plot the dephasing time Tss due to Vss 
for GaAs and Si as a function of the exchange splitting 
J. This dephasing mechanism is enhanced by larger J 
because of the J-dependence of the Vss coupling, which 
in turn reflects the fact that the Vss term originates from 
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FIG. 4: The timescale Tss of S-To coherence decay due to the 
Vss term for (a) GaAs DQD with Nl = Nr^10'^ at B = 200 
mT (l/ixl ~ 5.0 /ieV), and (b) DQD based on natural silicon 
with Nl = Nr = 10^ at B = 30 mT (I/^tI ~ 5.0 fj.eV). The 
diagonal hatching regions are where a± <C | J ± /it], |/iT| is 
not fulfilled. 



the J-induced asymmetry between the two virtual pro- 
cesses shown in Figure [1] Numerically, Tss is longer than 
Ta for smaller J, but shorter at larger J values, so that it 
can become the dominant dephasing channel for an S-Tq 
qubit in the large- J regime. We will discuss the crossover 
in more detail later in this Section. 

The suppression of Knight shifts experienced by the 
nuclei in a DQD (compared to a single spin) has impli- 
cations beyond simplifying the dephasing calculation in 
this subsection. For example, under certain conditions 
the transverse components of the Overhauser field can 
be treated by simply averaging the system's evolution 
over the distribution ()18p of classical Overhauser fields, 
with only dynamics from the Larmor precession of the 
nuclei (see Appendix [Ct . For a single spin in a QD, this 
approach is valid when the evolution time t is shorter 
than the inverse of the typical spread of nuclear spin 
Knight shifts,^^ i.e. t<7V^/yt.i^iSiZi Thus in this short- 
time limit the exact shape of the electron's wavefunc- 
tion is irrelevant. On the other hand, when t^ N,;, /Aa, 
this wavefunction shape has to be taken into account. 
The qualitative form of the decay then changes, and the 
semiclassical approximation (with the nuclei replaced by 
a classical field or a precessing large spin) fails. How- 
ever, for an S-Tq qubit in a DQD with a finite J and 
no 6't, the strongly reduced effective Knight shift for 
the nuclei means that the RDT calculations presented 
above are done in the "short time regime" for a single 
spin qubit, '^^'■-^ i.e. in the regime where the shape of the 
wavefunction is irrelevant. 

Another consequence of the suppressed Knight shift 
is the diminished importance of the bath dynamics in- 
duced by dipolar interactions among the nuclear spins. 
For a single spin qubit this dynamics leads to fluctu- 
ations of the qubit 's energy splitting by the longitudi- 
nal Overhauser field operator ^^Ailf. This dephas- 
ing mechanism causes a narrowed-state free induction 



11 



decay on a timescale of > 10 /is in a typical GaAs 
qq 30,36,72,75 (-^yjijig thesc predictions have not been yet 
verified, the theory given in these references correctly 
accounts for spin-echo decay measurements in GaAs^^ 
and sihcon^'''^^). For an S-Tq qubit with a finite J, the 
timescale of coherence decay due to dipolar dynamics of 
nuclei is expected to be longer by a factor of ^No^-t/A, 
which makes this dephasing mechanism irrelevant in the 
case of pure S-Tq dephasing. 



D. Dephasing due to Vse 

The Hamiltonian governing the V^e-induced dephasing 
takes the form H = Vse{\S) {S\ - \To) (To|). According 
to Eq. (I38|) . Vse consists of two independent (commut- 
ing) terms related to the two dots in the DQD, thus the 
total decoherence function is a product of the single-dot 
decoherence functions: WsB{t) = Wse,L{t) x WsQ^R{t). 

Vse contains both transverse and longitudinal compo- 
nents of the Overhauser field, so that the full RDT cal- 
culation is rather cumbersome (the semiclassical calcu- 
lation is also more involved than in the case of Vss, as 
shown in Appendix [C|) . On the other hand, it turns out 
that in most cases dephasing caused by this term is much 
slower than the previously considered mechanisms (due 
to Ha and Vss), and an exact calculation of coherence 
decay due to Vse is of little practical relevance. Below we 
present a calculation of a lower bound for the coherence 
time Ts0. This lower bound is obtained by replacing the 
59 operator in Eq. psp with (7^. The decoherence func- 
tion in this case is given by. 



Wse{t)^i:T, 



[pif 



(.^IoVso{r}dT 



T 



'^S^,Vse(T)dT 



(62) 

where T{T) is the time (anti-) ordering operator, and 
Vso^t) is in the interaction picture, with the nuclear spin 
operators given by Eq. (j5ip . 

Following a procedure analogous to the one described 
in the previous subsection, we obtain, for example. 



Tfe^(i) « Avga.^E^iAkAi e''^"'*^^ i 






t , (63) 



where the nuclear indices k and I refer to the nuclei in dot 
L and sinc(a;) = (sina;)/a;. The formula for T^ is analo- 
gous, albeit with a negative sign. The total decoherence 
function can then be obtained as 



simplified: 

g-iarctan(r;^j,t) ^i arctan(r;fj, t) 

Wse t«^-^M « . - , (65) 



where rj^g 



1 + te)" Vi + W) 



simplicity we assume a symmetric DQD (cr_L,L 
(Tj^/\/2), such that 



8wflCr^o-]^^y^. For 



O'.L.fl 



Wse U<w 



al3 



1 



1 + {mety 



(66) 



where i^se 
then given by 



^ve(7:,a\. 



The characteristic decay time is 



Tse> 



2J|mt| 



a^a . 



1 - ( J/a^t)^ 



2 - {J|^i^f 



(67) 



Although Eq. (p7|) is derived assuming a low i?-field 
in a hetero-nuclear material (so that T^b is shorter than 
w~„), it remains a good estimate of the characteristic 
decay timescale at higher fields in III-V materials. In 
Fig. [S] we show results of example calculations of Tse de- 
cay times for GaAs and Si. The non-monotonic behavior 
of Tse as a function of J is due to dependence of the cou- 
pling strength vg on J. As shown in Fig. [1] there are two 
terms contributing to Vse'- one associated with virtual 
transitions involving T_, the other involving T^. When 
J increases, the former is suppressed while the latter is 
enhanced, and a maximum in Tse (J) appears as a result 
of this competition. 
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FIG. 5: Lower bound for Tgg as a function of J for (a) GaAs 
with Nl ^ Nr = 10'' at B = 200 mT (/iT ~ 5.0 /leV) 
and (b) natural Si with Nl ^ Nr = 10*^ at B = 30 mT 
(I/itI ~ 3.5 fieV). The diagonal hatching regions are where 
cx <Si I J ± A(t|, |/^t| is not fulfilled. 



Nl 



-jarctan(A^ (t)) Nr ,; arctan(A-"(t)) 

Wse{t) = l[ I n / (64) 



™ ^l + {\l^{t)f n ^l + iX^{t)y 



,L(fl) 



L{B.) 



where Am are the eigenvalues of T^.^ matrices. 
When Wkit ^ 1, or for a homo- nuclear system, T^^ 



e,ki 



iyei^zy/o-kCiiAkAit, and the decoherence function can be 



E. Identification of the dominant dephasing 
mechanisms in GaAs and Si 

Now that we have investigated the individual dephas- 
ing channels for an S-Tq qubit by the hyperfine inter- 
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action, we would like to identify the dominant mecha- 
nisms in various regimes defined by the value of B field 
/iT , singlet-triplet splitting J, degree of nuclear reservoir 
narrowing np, dot sizes etc. Below we first list the ap- 
proximate formulas for the characteristic dephasing times 
due to each of the above-considered mechanisms. For the 
sake of clarity, in these formulas we assume Nl = Np — N, 
so that N]:) = N/2. Defining r = J/fiT, we have 



Tbc 
Ta 

Tss 
Tse 



4^6 



l + r 



a^ A 



> 



2Ve^- 

2J\tiT\ 



^l|i'-'-'i 



l-r^ 



(68) 
(69) 
(70) 
(71) 



Notice that while Tbc decreases with increasing J, even 
at r ~ 0.9 it is still much longer than N/A, which is « 10 
fis (~ 1 ms) in a typical GaAs (Si) dot with N « 10^ 
(10^). We thus do not include Tbc in the discussion 
below. 

Without nuclear state narrowing, there is a competi- 
tion between H^ and V55 induced dephasing. At small 
enough J, dephasing due to Ha dominates, while Tss 
becomes the shortest timescale at large J. We show ex- 
amples of calculations of all these times for GaAs and Si 
in Figs, m and [71 respectively. From Eqs. (|69l) and (|70|) 
we can obtain Jmax for which the two mechanisms give 
the same dephasing time: 



<-'lT 



0.64/iq 



(72) 



The approximate sign here can be replaced by an equal 
sign in the case of a single-isotope material (or when de- 
cay occurs at a timescale shorter than the Larmor pre- 
cession period of the nuclei). At Jmax the dephasing time 
is maximal: we anticipate a non-monotonic J depen- 
dence for the observed T2 time. The decoherence function 
should be well approximated by Eq. (|49)) for J <^ Jmax, 
and by Eqs. ^7^ and ^U^ for J> Jmax- 

When the nuclear state is narrowed with np < 1, Ta 
and (Tse) are extended by factors of n^ and nj, , respec- 
tively. With strong enough narrowing, the V^s-induced 
dephasing, i.e. dephasing due to hf-induced virtual tran- 
sitions between S and T± states, becomes the dominant 
source of hf dephasing in an S-Tq qubit. 



F. Dephasing in an InGaAs DQD 

We can apply the theory presented in this Section to 
two electrons in a vertically stacked self-assembled In- 
GaAs DQD, for which an estimate of the S-Tq coher- 
ence time has been experimentally determined recently42. 
Due to the strong interdot tunneling achievable in these 
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FIG. 6: (Color online) Dephasing times induced by Ha (solid 
red line) , Vss (dashed blue line) and Vse (dotted black line) as 
a function of J/ht for GaAs at (a) B = 1 T (|/jt| ~ 25 /ieV) 
and (b) B = 200 mT (I/xtI « 5.0 /ieV), both with Nl = 
Nb_ — 10^, resulting in az — o-± ~ 0.1 /leV with no narrowing 
{np^l). 
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FIG. 7: (Color online) Dephasing times induced by Ha (solid 
red line), Vss (dashed blue line) and Vse (dotted black line) 
as a function of J/ht for Si at _B = 30 mT (I/^tI ~ 3.5 /^eV), 
both with Nl ~ Nr = 10^, resulting in az = u± ^ 1.5 neV. 



structures, the values of J are much larger than those 
in gated GaAs DQDs. Most interestingly, one can in- 
vestigate the S-Tq coherence at the value of interdot de- 
tuning where J is to first order insensitive to the fluc- 
tuations of the electric bias (i.e. the charge noise), while 
still satisfying J ^ az- In Ref. \^ at such an "opti- 
mal point" with respect to charge noise, Jop « 100 /xeV, 
which is in fact larger even than the /it splitting due 
to the applied B field. While we have confined our- 
selves so far to the regime where J<\^t\, all our results 
are applicable for J > j/ixl (as long as we avoid the re- 
gion of strong hyperfine-induced S-T± mixing - in other 
words I J ± /ixl ^ CTj^ has to be satisfied). In Fig. [5] we 
present the decoherence functions WA{t), Wss{t), and 
Wse{t) for a Ino.5Gao.5As DQD. The dot parameters are 
Nl^Nr = 10^, J=100 /ieV, and 5 = 200 mT (the elec- 
tron g-factor is approximated by ^cff ~ 0.5). For these 
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parameters Ha and Vss are of similar importance, and 
Vss dominates the long-time decay. The calculated co- 
herence time is of the order of a microsecond, which is 
in qualitativ e ag reement with the lower bound of 200 ns 
given in Ref. |43J. 
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Hamiltonian in the new eigen-basis. Nevertheless, the 
physical picture is quite clear here. The mixing of S and 
To states means that in the new eigenstates, electron spin 
density in each dot does not vanish anymore. As such 
the linear longitudinal Overhauser field 69 would lead 
to dephasing between the eigenstates, similar to what 
happens to single spin qubits. Indeed, if ^t ^ J, the 
eigenstates would approach the product states again, so 
that we should recover dephasing of two independent 
spins exactly. In this Section we focus on this transi- 
tion from reduced dephasing in the highly-symmetric S- 
Tq two-level system to the single-spin dephasing in the 
large-field-gradient limit. 

A. S-To free-induction decay dephasing in the 
presence of an interdot magnetic field gradient 

At finite 9t, the eigenstates of the first matrix in 
Eq. (dni) are 

1^') = cos7|5)+sin7|ro) 
|r^) = -sin7|5)+cos7|ro). 



FIG. 8: The decay of the decoherence function, calculated 
separately for each of the mechanisms, for Ino.5Gao.5As DQD 
with two electrons at 5 = 200 mT (assuming the g- factors of 
the electrons in both dots to be «0.5), with J =100 /leV and 
for Nl^Nr^ 10^ Here \ht\ r; 5.8 /ieV. 



IV. CALCULATIONS OF S-To DEPHASING IN 

THE PRESENCE OF AN INTERDOT 

MAGNETIC FIELD GRADIENT 

In the presence of a finite field gradient 9t ^ o^z (ei- 
ther from a nanomagnet, or due to a previously prepared 
Overhauser field gradient), one needs to obtain the new 
eigenstates that account for the 0T-induced mixing of S 
and To states, and then to re-derive the pure dephasing 



where the mixing angle 7 is defined by 
tan 27 



29t 



(73) 



The full range of 7 is from to 7r/4, corresponding to 
limits oi 9^ <^ J and 9^ ^ J, respectively. The corre- 
sponding eigenenergies are 
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(74) 
(75) 



The complete dephasing Hamiltonian in the basis of 
5' - T^ is then 
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,56'sin27 
59 cos 27 



59 cos 27 
-(50 sin 27 



Vss cos^ 7 -I- Vh sin 27 



-Vss 



sin 27 



Vu cos 27 + Vah 



-Vss 



sin 27 



r 



Vh cos 27 — Vah Vss sin 7 — Vh sin 27 



(76) 



where we have neglected Hb and He since their influence 
on dephasing is very weak, as discussed before. 

The influence of the off-diagonal terms in Hgrp could 
be included by performing another canonical transfor- 
mation. However, the off-diagonal terms related to the 
transverse Overhauser fields generally only give sublead- 
ing corrections to the diagonal terms, except at the 



9t ^ J (when the diagonal Vss term vanishes) or 
9t -^ J (when the diagonal 59 and Vh terms vanish) 
limits. Keeping only the lowest order terms (up to sec- 
ond order in hf coupling strength), we obtain the pure 
dephasing Hamiltonian in the S' — Tq space as 
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H 



STo 




sin 27- 2^008^27 






'sin 27 + ^008^27 



Vss cos^ 7 + Vh sin 27 ' 

Vss sin^ 7 — Vh sin 27 y 



(77) 



r 



We can define H'^ — —59'^ cos"^ 27/ J, which goes back to 
Ha when 6't = 0. 

The key new feature here is the re-appearance of the 
(50-hnear terms in the dephasing Hamihonian. The 59- 
induced inhomogeneous broadening that plagues single- 
spin qubits is now back in action, though its effect is 
reduced when the field gradient 9^ is small compared 
to the exchange splitting J (when the mixing angle 7 is 
small). The inhomogeneous broadening dephasing time 
due to the ^0-linear terms is given by 

V2_ V2J ^,3) 






1 



4CT.07 



I sin 27 1 (72 

where the approximate formula holds when 7 <C 1 (but 
for 9t still larger than az)- When 7 approaches 7r/4 
(i.e. 9t^J), 2^2* St approaches the Tg ~ I/ctz for a single 
spin in a QD, which is about 10 ns for a typical GaAs 
QD [see Fig. ^]. On the other hand, for 9t < J, T*^g^ ex 
J/9tCz- As 9t decreases, Tj g grows, until we reach the 
regime 9^ ^ CTj, when dephasing due to the quadratic 
H'j^ term becomes more important. In general H'^ leads 
to a characteristic dephasing time of 

„2 



T\ = 



J 



a^ cos-' 27 



(79) 
1, and 



and T^ due to the 



When 0T ~ CTz ^ J, 7 ^ 1, so that cos 27 
T'^K. Ta- In Fig. IHwe plot T* 
longitudinal Overhauser fields in response to varying ^x 
at fixed values of J . A transition of T^ g from two-spin 
to single-spin coherence dynamics is shown in Fig. ^a) 
as 0T approaches J/ 2. When 9^ ^ cr^, the resulting 
two-spin dephasing is essentially the same as in the pre- 
viously discussed case of 9t = 0. Instead of T2 ,T^^Ta 
now represents the longitudinal Overhauser field induced 
inhomogeneous broadening. 

The remaining terms on the diagonal are related to the 
previously discussed Vss and Vsg terms. The dephasing 
due to the first of them 



Vss cos^ 7 





Vss sin^ 7 



(80) 



is calculated in the same way as in Section llll CI (compare 
also with the Hahn echo calculation in the next subsec- 
tion), only with T^f} term in Eq. (j56p multiplied by cos 27. 
The resulting dephasing time is 

Tss 



T' 

^SS 



cos 27 



(81) 
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FIG. 9: The dephasing time (a) 72*9^ due to the (5^-linear 
terms and (b) T'j^ from the second canonical transformation 
in GaAs as a function of 9t at two fixed values of J. The rms 
of the Overhauser field difference between the dots is CTz =0.1 
/.ieV. For 6t 2> J, T2*,9t ^^ saturated at ~ 10 ns. 



which shows that the dephasing due to the transverse 
Overhauser fields described by V55 is suppressed by the 
presence of the field gradient 9^. On the other hand, the 
influence of the interactions previously appearing in Vse 
term is enhanced by the finite 0t , similar to the inhomo- 
geneous broadening due to 59. The dephasing due to the 
term 



Vh sin 27 

-VHsin27 



(82) 



is calculated in a way analogous to the one in Section 
IIIIDi only with az in Eq. ([63| replaced by 9^. The co- 
herence dynamics is thus described by Eqs. (p5)) and (p5)) 
in which rjss is replaced by ijh = —'riseJ siTi.2j/2az- The 
resulting dephasing time is given by 



Th^ 



'2(7z Tse 
~|sin27| 



(83) 



When 9t <^ J , Th ^ §^Ts9, with Tse given by Eq. §7}. 
Clearly, this dephasing channel is strongly enhanced 
when 0T ^Cz. As 9t increases, so that 9t ^ J ^ ctz, 
Th -^ 2Tse(Tz/J ^Tse- In Fig. [10] we show examples of 
calculated values of T'gg and Th- One can see that with 
exception of the regime of small 9t, these decay times 
due to the transverse Overhauser fields are much shorter 
than the T^ g time from Eq. ([75| shown in Figure IH^. 
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FIG. 10: The dephasing times Tgg from Eq. dHJ and Th 
from Eq. ^^, due to the terms given in Eq. ((MJ and Eq. ^^ 
respectively for GaAs DQD with Nl^Nr^ 10** at B = 200 
mT {\iit\ Ri5.0 peV). 



B. Spin echo for the S-Tq qubit in the presence of 
an effective magnetic field gradient 

A finite field gradient (i.e. Ot > (Jz) can produce con- 
trolled rotations around the x axis of the S-Tq Bloch 
sphere when J = 0. Such rotations allow faster prepa- 
ration of superposition of S and To states (see Appendix 
IX)) . but more importantly, they allow the Hahn-echo 
experiment^ in which a 7r-rotation about the x axis 
[from here on we will label a rotation of angle around 
the a-axis as a ((?!>)a pulse] is performed in the middle of 
a free evolution (i.e. at T = t/2). If a qubit is initialized 
in the S state, the full Hahn echo sequence would con- 
sist of a (7r/2)a;-pulse (to prepare a S-Tq superposition), 
free evolution for time t/2 in the finite- J regime, then a 
TTx pulse, another free evolution period of t/2, and finally 
another (7r/2)2:-pulse, which prepares the qubit for the fi- 
nal projective measurement on the S state. After simple 
manipulations (similar to the ones shown in Appendix 
E]) we arrive at the singlet probability at the end of the 
sequence: 
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(84) 

where H is the total Hamiltonian in the S-Tq subspace. 
Now we rotate to the S" — Tg basis. As we did previously, 
we again make a pure dephasing approximation for the 
Hamiltonian H' in the rotated basis, such that 
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(85) 
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is the truly "echoed" part of the signal, while 
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(2r) = Tr/ (p/e 
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(86) 



(87) 



is the component of the signal which decays just like 
the free evolution signal. The presence of this term is 
due to the fact that the field gradient used for rotations 
about the x axis is "always on" in the current-generation 
experiments 42 Here H'g and H^ are the diagonal terms 
of the S' — Tq pure dephasing Hamiltonian. Using a ro- 
tated coordinate system on the Bloch sphere (with -|-z' 
direction corresponding to a rotated 5" state), all the 
rotations are slightly tilted away from x' axis. The re- 
sulting pulse error leads to an appearance of a part of 
the signal which corresponds to evolution unaffected by 
the TT pulse. 

Note that for \29t/J\ < 1, 27 < 1, and the FID- 
like part of the signal is just a small correction. Only 
a.t 2 9t ^ J the two components contributing to Ps (t) 
would be comparable. Since it is expected that W^^{t) 
exhibits slower decay than W^^^{t), in this regime the 
echo decay time will be close to the FID decay time cal- 
culated previously. Below we focus on the J » 2\9t\ 
regime, and calculate the time dependence of W^^{t). 



Echo decay due to longitudinal Overhauser fields 



The Hahn echo sequence cancels completely any con- 
tribution to dephasing of H'^ q terms. Furthermore, the 
phase accumulated during the evolution due to the diag- 
onal 59 sin 27 and if^ terms is completely cancelled. 



2. Echo decay due to the transverse Overhauser fields 

There are two terms in H' due to transverse Over- 
hauser fields. The first is the modified Vss interaction 
term: 
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(88) 



where /' and ct^ are identity and Pauli-z matrix in the 
S'—Tq subspace. The identity part of the interaction does 
not cause decoherence between the S' and Tq states, and 
will be dropped from the following considerations. The 
second is the large-0T analogue of Vse term considered 
previously: 



V^ = Vu sin 270-!, 



201 



Vnij' 



(89) 
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Below we will consider their independent contributions 
to dephasing. 

The calculations in both cases are straightforward 
modifications of spin echo theory of Refs. I3lll32ll34l . Fol- 
lowing the derivation of Ref. |32|, the W^^{t) function 
from Eq. (|86p can be written as 



W"^{t)=TTi 



piTce^p{-t lj{QV{QAt'c 



(90) 

where Tc is the operator ordering V'{tc) on the closed 
time-loop contour, -^^ ^'(^c) is the respective interaction 
term [from Eq. dMl or Eq. ((Ml)], tc = {t,c) is the time 
variable on the contour (with c = ±), and /(ic) is the 
temporal Hahn echo filter function defined on the con- 
tour. 

As discussed previously, Vgg{t'^) interaction is given by 
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i 27 Vss ^ < 
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where 0i = 1 for i (£ L and (fi 
corresponding filter function is 
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-1 for i e R. The 



/(t^) = e(t - t')Q{t' - 1/2) - e(t/2 - t')e{t') . (92) 

where 9(x) is the Hcaviside function. The resulting T 
matrix is 
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where No = l/(iV^^ + iV^^). Following Ref.lH, leads to 

^ (94) 



Wg'^it) 



l + \Rss{t) 



with 



AlAl sm 



4 ^apt 

4 



Nl 



'^ap 



Rss{t) = 16wL cos^ 27 ^ aaapnaUp 

(95) 
We thus see that the Hahn echo S-T^ coherence decay 
due to Vss interaction is of the same form as the echo 
decay of a single spin coherence (cf. Ref. ISlUsa ). 
For Vg{t'^) interaction we have 



Veit'c 



k,l 



. (96) 
and the filter function is given by Eq. (|M|) . Since V'g is a 
sum of two commuting terms (under our approximation 
of neglect of the overlap of L and R wavefunctions even 
at finite J), W^^{t) is a product of two functions, each 
corresponding to one of the dots. Assuming symmetric 
dots, i.e. Nl—Nr = 2Nd we have then 
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(97) 
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FIG. 11: (a) Wss (t) and (b) Wg'^{t) for GaAs DQD with 



Nl=Nr = W'' &t B = 200 mT and i 
values of J. Here I/itI ~ 5.0 /^eV. 
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FIG. 12: (a) Wgs {f) and (b) Wg''{t) for GaAs DQD with 
Nl ^ Nr = 10® at B = 700 mT and St = 2.5 mT with various 
values of J. The parameters here are chosen to correspond 
closely to the ones employed in recent experimental work from 
Ref. \M- 
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(98) 
The decoherence functions due to both of the above in- 
teractions are shown for GaAs in Figures [TT] and [T^l It 
is obvious that the Hahn echo signal shows more pro- 
nounced oscillations at lower B fields fFig. fTTj) . when the 
processes of virtual fiip-flops between the 5- To subspace 
and the polarized triplets have more importance. Let us 
mention that in a recent experiment^ on Hahn echo in 
the large J regime, a decay of echo signal on timescale of 
a few microseconds was seen at B = 0.7 T. Comparison 
with results presented in Fig. I12[ which show only visi- 
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bility loss of less than one percent even at very large J, 
clearly supports the main claim of Ref . \^ that the charge 
noise (specifically the fluctuations of J) is the dominant 
source of dephasing in this experiment. 



V. EXCHANGE GATE ERROR 

So far we have focused on hf-induced decoherence for 
an S-Tq qubit with a finite exchange splitting J. In the 
case of single-spin qubits, exchange interaction plays the 
crucial role of generating two-qubit gates such as SWAP 
and ControUed-NOT gates.— In this Section we evaluate 
how hf-interaction affects the fidelity of such two-qubit 
gates. To allow a simple SWAP gate between two single- 
spin qubits, we assume there is no field gradient between 
the double dot: ^t = 0. 

Let us first make a qualitative examination of the gate 
error problem. Recall that the magnitude of the exchange 
splitting J is generally between 0.1 /zeV and 10 ^eV, 
corresponding to a gate time Tg between 40 ns and 0.4 
ns. Comparing the short gate time with the decoherence 
times we have calculated in the previous two Sections, the 
only factor that could significantly impact the fidelity of 
a SWAP gate within the S-Tq subspace is the fiuctuations 
in the Overhauser field difference, S6. In the absence of a 
field gradient, S-Tq states are affected by 69 only at the 
second order level through Ha , as we have shown in Sec- 
tion [1111 On the other hand, the dephasing Hamiltonians 
for the polarized triplet states, Ht+ and Ht_ (given by 
Eqs. (|E3|) and JES} in Appendix |E1 respectively), do con- 
tain (5/i, the fluctuations in the average Overhauser field. 
One can reason then that the SWAP gate fidelity would 
suffer much more severely if polarized triplet states are 
involved. 

We focus on the fidelity of a SWAP gate in a two- 
electron DQD, which is defined as 



F = (l/'in|t^sWAP Pout t^sWAplV'in) 



(99) 



where pout refers to the output density matrix and C/swap 
denotes the SWAP operation. The two qubits are initial- 
ized in JT/iin) = iV'i) (8) \ip2), where |V'j) = a^ It) + h \i) 
with i = 1,2. Note that the gate fidelity F defined here 
is input-state-dependent. Furthermore, it can be easily 
generalized to the cases where the input state is repre- 
sented by a density matrix instead of a pure state. In the 
current study we only investigate a few representative 
examples instead of giving a comprehensive discussion. 
However, the results we present should already paint a 
clear picture on how hyperfine interaction affects the ex- 
change gate fidelity in general. 

Here we focus on three types of input states: (1) a fully 
polarized state, for example jV'in) = |tt); (2) an initial 
state in the S-Tq subspace, for example \tjjin) = |ti); and 
(3) a general initial state that contains both the polarized 
and unpolarized two-spin eigenstates. 

We start with a fully-polarized pure input state, 
l^in) = \T+) or |T_). For either of them, the gate fi- 



delity always stays at 1, since H^f does not mix two spin- 
up (spin-down) states, and the exchange gate does not 
change the state at all. If the input is a superposition of 
the two polarized triplet states, such as (|T-|_)-|-|T_))/-\/2, 
the exchange again gate does not change either of the 
triplet states, though the two states do dephase relative 
to each other because of the S^ term in Ht± , irrespective 
of whether the exchange gate is turned on or off. 

When we have an input state in the S-Tq subspace, 
it can be expressed as a linear combination of \S) and 
|To) with appropriate coefficients. For |^in) = \ti) = 
{\S) + \To))/V2, the gate fidelity Fi takes the form 



Fi 



1 1 



ReWsToit) 
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which also applies to the case of |?/'in) = lit)- We have 
studied the decoherence function Wstq (t) in detail in Sec- 
tion |TIT1 which decays due to the quasi-static longitudi- 
nal Overhauser field (Ha) and due terms proportional to 
transverse Overhauser fields, Vss and Vse)- When the 
nuclear bath is not narrowed. Ha = —0^/J is the lead- 
ing dephasing mechanism for smaller J (cf. Figs. [6] and 
[7] for GaAs and Si, respectively). The resulting gate fi- 
delity can be approximated as (for simplicity, we assume 
a square pulse of J with a SWAP gate time of Tg — it/ J) 



Fi 



1 1. 
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2 2 
and the gate error is 
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(102) 



There are two interesting features here. First, the gate 
error is inversely proportional to J*. Thus the larger 
the J, the smaller the gate error (although the charge 
noise induced error would more likely be the dominant 
factor for larger J values, and J should be smaller than 
"/max ~ 0.64^T in order for iJ^-induced dephasing to 
dominate over t^s^-induced one, see Section |IIIE|) . Sec- 
ond, narrowing of the nuclear bath (i.e. reducing np) can 
help the gate fidelity quite dramatically 1 — Fi ex rip . In 
Fig. [13] we plot the SWAP gate error 1 — Fi with the input 
state IV'in) = Itt) as a function of the exchange splitting 
J. Both the strong J dependence and the strong np 
dependence are clearly illustrated. 

In more general cases, IV'in) contains both polarized 
and unpolarized two-spin states (and it does not have 
to be a product state either). Dephasing between the 
polarized triplet states and the unpolarized states be- 
comes relevant. Different from the Sz = subspace, 
where Ha is the leading cause for dephasing, here it is 
the average Overhauser field Sfi that leads to inhomo- 
geneous broadening and dominates dephasing. In GaAs 
this inhomogeneous broadening leads to a Tj* time in the 
order of 10 ns, which trumps all other dephasing mech- 
anisms due to higher-order terms in the effective Hamil- 
tonian from Eqs (|E2IIE5|) . For example, for an initial 
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state IV'in) = Tjlt + D^-Tflt — i); using the classical 
avergaing procedure we employed in previous Sections we 
obtain the SWAP gate fidelity as 



Fo 






1 






(103) 



where T^ is the characteristic time due to the inhomoge- 
neous broadening effect from 6(1, 



T,= 



V2 



V2 



i<L + <r) "^^^ 



(104) 



Here we have assumed that a^ — \ /cr? 



i.e. we include the possibility of the narrowing of the 
Overhauser field distribution in each of the dots sepa- 
rately. Note that this is distinct from the previously con- 
sidered case of narrowing of Overhauser field difference 
between the two dots, leading to our introduction of a 
narrowing factor n'p instead of np used above. Narrow- 
ing of this difference (i.e. decreasing (Jz, the standard 
deviation of the field difference, from its natural value 
by a factor of np) can occur when no narrowing of the 
field distribution in each of the dots is present - only a 
correlation between the z components of the Overhauser 
field in L and R dot is needed. Keeping this in mind, we 
see that the gate error takes a simple form of 



T2 
1 _ ^2 - — 



7r^(n 



Fl 



J2 



(105) 



Since a± is in the order of 0.1 jieV for GaAs, the ex- 
change coupling J needs to be larger than 1 /xeV for the 
gate error to be manageable in the context of quantum 
information processing. In Fig.[T3]we show our results of 
1 — F2 as a function of J for the input state given above. 
Comparing to the curves for the initial state in the S- 
Tq subspace, the gate error here has a larger magnitude 
and a weaker J-dependence. This difference comes from 
the fact that 1 — Fi is dominated by Ha = 69^/ J, while 
1 — i^2 is dominated by Sfi in Ht± . 

In short, while high-fidelity exchange gate can be ob- 
tained when an input state is in the S-Tq subspace, gate 
error is generally larger for an input state containing po- 
larized states, dominated by the inhomogeneous broad- 
ening from the quasi-static nuclear spin reservoir. Since 
the later type of input states is the norm rather than the 
exception, the upper two curves in Fig. [13] are more rep- 
resentative of the qualitative behavior of a general input 
state. Furthermore, the numerical coefficient for the gate 
error 1 — i^ is input-state-dependent: it is determined 
by the weight of \T±) within the input state. Finally, 
it is also important to point out that our consideration 
in this Section assumes the simplest implementation of 
the SWAP gate. As pointed out in Refs. ItsIItqI . pulse 
shaping and other means can be employed to reduce the 
sensitivity of an exchange gate to the Overhauser field 
fiuctuations. 
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FIG. 13: The SWAP gate error in GaAs at B = 200 mT 
for two input states: (a) l^in) = IV') ~ Iti) s^nd (b) ji/'in) = 
\^) = -1= (It) + |;» ® ^ (It) - ID). Solid curves show the 
gate errors with the Overhauser fluctuations a^ ~ 0.1 /xeV, 
while dashed curves are with az « 0.05 fieV. Here az — n'pa±, 
as discussed in the main text below Eq. (|104p . 



VI. SUMMARY AND CONCLUSIONS 

With the aim of helping to understand experiments 
related to S-Tq qubits and single-spin qubits, we have 
performed a theoretical study of hyperfine-interaction- 
induced dephasing for two electron spins in a dou- 
ble quantum dot at finite S-To splitting J. For the 
S-Tq qubit, we evaluate various hf-induced dephasing 
terms between the singlet S and the unpolarized triplet 
state To, and identify the dominant dephasing channels. 
Specifically, we find that in the absence of a magnetic 
field gradient across the DQD, inhomogeneous broaden- 
ing for an S-Tq qubit is significantly suppressed relative 
to a single spin qubit, with the most significant source 
of inhomogeneous broadening being the Ha term, pro- 
portional to the square of the fluctuations of the longitu- 
dinal Overhauser field difference between the dots. An- 
other important decoherence channel is the Vss term, 
which originates from dressing of the singlet state by 
the polarized triplets due to transverse Overhauser field. 
Throughout a wide range of parameters, these two terms 
compete for the role of the dominant source of decoher- 
ence, leading to a non-monotonic dependence of the T2 
dephasing time on J. 

The physical picture changes when there is a finite 
magnetic field gradient across the double dot. The gra- 
dient dictates that the two-spin eigenstates are now su- 
perpositions of the singlet and unpolarized triplet states, 
and there is a finite electron spin polarization in each 
of the two quantum dots. Random longitudinal Over- 
hauser field difference can now cause dephasing between 
the two-spin eigenstates, in analogy to inhomogeneous 
broadening in a single spin qubit. In other words, the 
S-Tq encoding acquires more of a single-spin-qubit de- 
coherence characteristics as the singlet-triplet mixing in- 
creases. This increased hf-induced dephasing is the price 
that one has to pay for having two-axis control (with the 
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X axis rotations provided by the field gradient) over the 
S-Tq qubit. 

We have also considered the Hahn echo experiment 
in which the superposition of S and Tq states evolves 
at finite J, while being subject by a 7r-pulse (generated 
by a built-in gradient of the effective magnetic field be- 
tween the dots, which drives the rotation when J is sup- 
pressed for a chosen time) in the middle of the evolution 
period (for a recent experiment see Ref. |42| ). The dy- 
namics of the echo signal is then due to the transverse 
Overhauser field terms, which lead to appearance of os- 
cillations (closely related to the ones predicted^i^^ and 
observed^^ in the J w regime) , the amplitude of which 
increases as the applied magnetic field is lowered. The 
observation of such oscillations in a Hahn echo experi- 
ment in the large J regime would be a clear signal of 
achieving a strong suppression of the charge noise (caus- 
ing fluctuations of J) in the system. A non-monotonic 
dependence of the characterstic free-evolution decoher- 
ence time T2 on J (with the maximal coherence time seen 
at Jmax~0.64/iT) is another predicted experimental fea- 
ture signifying the dominance of hyperfine mechanisms 
of decoherence in the system. 

We have also studied hf-induced gate errors in an 
exchange-gate for two single-spin qubits, with a particu- 
lar focus on the SWAP gate. Due to the short durations 
of an exchange gate, inhomogeneous broadening from the 
longitudinal Overhauser field is the dominant source of 
gate errors. By exploring several typical initial states, 
we are able to obtain the dependence of the gate errors 
on J and the standard deviation of the distribution of 
the Overhauser fields in the two dots. The obtained re- 
sults show that for realistic multi-qubit operations with 
single-spin qubits, the fluctuations of longitudinal Over- 
hauser fleld remain a serious source of computation errors 
and control schemes focused on mitigating these errors by 
exploting the quasi-static nature of the Overhauser field 
fluctuation, similar to the ones presented in Refs. I78ll79l 
for single qubit operations, should be seriously consid- 
ered. This necessitates a generalization of the currently 
proposed techniques for noise-resistant single qubit oper- 
ations to the multi-qubit case which is still very much in 
its infancy)^ 

While our calculations are focused on two coupled spin 
qubits, the lessons we have learned should be useful for 
analysis of decoherence of multiple coupled qubits. For 
example, our results on the S-Tq qubit in the absence of 
fleld gradient indicate that electron interaction and sym- 
metry can help suppress decoherence, in the same man- 
ner as decoherence-free subspace, even though contact 
hyperflne coupling is a completely local interaction. On 
the other hand, if the symmetries are broken, whether 
by intrinsic or extrinsic inhomogeneities, the coherence 
properties of the overall interacting system is more simi- 
lar to an ensemble of individual components. Neverthe- 
less, based on our results for S-Tq decoherence, we may 
speculate that dephasing between two states of a cou- 
pled n-qubit system should generally be slower than two 



states of an uncoupled n-qubit system. 
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Appendix A: Possible w^ays to measure the S-Tq 

coherence at large exchange splitting J in gated 

DQDs 

With the full control over the S-Tq qubit available;^ a 
measurement of coherence between S and Tq states can 
be performed at nonzero J in the following way (which 
in fact has been recently realized in Ref. |42) . Using the x 
rotations (enabled by a nanomagnet located close to the 
DQD^— i^i^ or a pre-established gradient of the nuclear 
polarization between the two dot o^^'^^ ) one can rotate 
the S state initialized in a standard way^ into a super- 
position of S and To (such as |±r) = {\S) ± i |To))/V^ 
states), and then set J to be much larger than the 6't 
term. At such large J the infiuence of the environment 
is expected to lead only to pure dephasing of the super- 
position of states - the diagonal elements p^^ and p^ ^ 
of the reduced density matrix of the qubit are conserved, 
while the off-diagonal elements, p^^^ = Pt s ^^^ decay 
due to fluctuations of the energy splitting of the S and 
To states. After letting the state freely evolve for time t, 
one can then measure the degree of coherence (i.e. p^^ ) 
by performing the state rotation which would bring the 
qubit back to the S state provided that the actual state 
after time t is the same as the initial one. Subsequent 
measurement of the singlet return probability P5 (i) is a 
measure of how much dephasing occurred in time t. 

The coherence-measurement protocol can be formal- 
ized as follows, assuming a constant finite 6't. The qubit 
is first initialized into the singlet S state<^ At t = 0, 
a {■7t/2)x pulse is applied (by pulsing the DQD to the 
J <ti 9^ regime and stay there for a period of time r, 
such that 6'tt = 7r/4), so that the qubit is rotated to 
the initial state of | — y). The DQD is then kept in the 
large- J regime and evolve freely for a time period t. Af- 
ter this time another (37r/2)j, rotation is performed, and 
subsequently the singlet probability Ps{t) is read out by 
charge sensing using the Pauli spin blockade. I Assuming 
that all the gate operations here take negligible time com- 
pared to the free-evolution time t, the singlet probability 
Ps{t) is defined as 



Ps{t)=TT[\S){S\m] , 



(Al) 
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where the final density matrix is given by 



m = U: 



37r/2,2;f 



~iHt 



Pq{Q) ® pi{0)e'"'U^ 



37r/2,: 



(A2) 



with /5q(0) and yO/(0) being the initial density matri- 
ces of the qubit and the nuclear bath, respectively, and 
Ux,3tt/2 = — (l+^o'a:)/ v2 being the unitary transformation 
corresponding to the final rotation of the qubit 's state. 
After simple manipulations we obtain 



Ps{t) 



1. 



rTr, 



(/5/(0) 



e'"'a,e''"'a. 



(A3) 



Since the DQD Hamiltonian is a pure dephasing one at 
large J,H^Hs \S) {S\ + Ht, |To) (To| (where Hs, Hto 
are operators in the Hilbert space of the nuclear bath) , 
we obtain 



Ps{t) = l + lReWsToit), 



where 



WsToit) 



_ Psnit) 



P^stM 



Tr/ (p/(0)e'^^o 



t^-iHst 



(A4) 



(A5) 



In other words, the time evolution of the singlet proba- 
bility Pg (i) carries the complete information of the S-Tq 
coherence Wstq {t) ■ 

In the absence of a controlled 6t gradient field, it is still 
possible to create an initial state being a superposition of 
\S) and |To) by intializing the S'(0, 2) state, and then adi- 
abatically lowering J, so that the state obtained at J <^az 
in the (1,1) charge regime is an eigenstate of the electron- 
nuclear-spin Hamiltonian. Depending on the sign of the 
expectation value of 69, one of the \±X){\S) ± \To))/V2 
states is created. The S-Tq splitting can then be rapidly 
increased to a finite value J, and the DQD system is al- 
lowed to evolve freely for time t. After this evolution J 
is again rapidly reduced to very low values, such that the 
state is projected onto the two-electron product state ba- 
sis. Subsequently the DQD is adiabatically swept to the 
spin blockade regime for the standard readout procedure 
via spin blockade. A calculation similar to the one given 
above gives us 



1 



^s(i) - o + TTr/ /5/(0)a,e 



CTa-e 



(A6) 



If we now make the pure-dephasing approximation for 
the Hamiltonian, we arrive again at Eq. (|A4p . 



Appendix B: Derivation of hyperfine coupling 
Hamiltonian 



To focus on the coupled electron-nuclear-spin dynam- 
ics, we project the hyperfine interaction on the lowest- 
energy two-electron states and obtain Eq. ((4]) in the main 
text. The basis states are the ground singlet and triplet 
states: 



I c,\ in) - i;t) ^ , , , 

\S) = 1= (E) -0s(ri, r2) , 



\To) = 



V2 

in) + i;t) 

V2 



'^As(ri,r2): 



\T+) = \n)<E>iJAs{ri,r2), 
|T_) = |;;>®^As(ri,r2). 



(Bl) 



Here '0(A)s(ri, 1*2) denotes the (anti-) symmetric combi- 
nation of the orbital envelopes ^^(r) and 'i'ij(r). 



V's(ri,r2) 



^As(ri,r2) 



[*L(ri)*fl,(r2) + *i?(ri)*L(r2)] 



i^EVW+W) 



(B2) 



[*L(ri)*fl(r2) - *7?(ri)*L(r2)] 



W2(1-|XP) 



(B3) 



with the orbital overlap x = ir /^'^^K'^)*^!'") ' '^^^^ 
1^0 being the volume of the primitive unit cell (note that 
we use the normalization in which the envelope wavefunc- 
tions are dimensionless, and all the dependence of the hf 
interaction on the periodic parts of the Bloch functions is 
carried by Aa constants below). Generally x <C 1, since 
the local orbitals decay rapidly away from the center of 
each dot. 

To calculate the hf matrix element, we need to first 
evaluate the matrix elements of the 5-functions. Specifi- 
cally, 
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{'4's\Aai{\S{rk -ROIV^As) 

Aa[{\voS{rk -Rj) 



= I dri dr2 



[*I(i-i)*fl(r2) + ^fl(ri)*i(r2)] X [^L(ri)*H(r2)-*fl(ri)*L(r2) 






{|vI/i(ROP - |*fl.(R.)l'-2Im[x*n(R.)*i^(R.)]} ^ (-l)'^'+^i?.(x) ,fc = 1,2 . (B4) 

1 



2(1 -IxP) 



|vI/i(R,)|2 + |vI/j,,(R,)|2_2Re [x**I(R.*fl(R.)]} (B5) 



(^sl A^i^i^oSivk ~ R.) I^s) = A(x) = a(x) + 2(1—^^^ [x**HR^)*fl(R. 

I 



(B6) 



In contrast to Bi{x), C'i(x) and Di(x) are both symmet- 
ric in (ri,r2). 

Writing Sfe • I, = S-^Jf + i(5+/,r + 5'^:/+), we can now 
evaluate the hf matrix elements straightforwardly. For 
example, 



{S\Hm\To) = J2 {^JslJ^M^^o^i^k 



RO 1^. 



'as; 



fe=1.2 



., (ni s^i^ in) - gti s^jf i;t) 

2 



The generalized form of the complete hf coupling matrix 
is then: 



^hf(x) =, 



/ 

B^ix)^ 



V2 

Mx) 

V2 



B,ix) j+ 



B^{X)^ 



cAxIt- 

V2 ' 

g.(x) 7-+ 

V2 i 



BiM/+ EiMjr\ 



V2 

g.(x) 7-+ 

72 ^i 

Q(x)/f 





V2 

c'dx) 




c.(x)/f/ 

(B7) 

In Eq. (jB7l) , we can neglect terms containing x and/or 
^i(Ri)^|j,(Ri) . This approximation is justified because 
the number of nuclei in the overlapping region is very 
small, and their interaction with the electron very weak 
(as they are located in the tail regions of both orbitals). 
After some simple algebra we find, for instance. 



{S\ Hm ITo 



EjfA' 



A?)l: =Y.B,Ii - (B8) 



The rest of the hf matrix elements can be obtained in 
a similar manner. With a negligible x i the hf coupling 
matrix can be approximated as Eq. ^. 



Appendix C: Semiclassical approach to the 

Overhauser fields and the dephasing induced by 

averaging over them 

The slow dynamics and the large numbers of nuclear 
spins in a quantum dot makes a semiclassical description 



of their average effect an enticing prospect. Indeed, it 
has been pointed out'^* that at the short-time limit, be- 
fore the electron-mediated nuclear dynamics causes any 
significant effects on the nuclear spin dynamics, a semi- 
classical description that account for only the Larmor 
precession of the nuclear spins is sufficient to explain 
most experimental observations. Here we give a brief 
summary of the semiclassical picture of the Overhauser 
fields, and estimate the magnitude of various hf-related 
terms appearing in our effective Hamiltonians. 

We write the Overhauser operator as h.L/R — 

'^- A^ 3i and treat its mean field average as a clas- 
sical field h^//j. The transverse Overhauser operator is 



given by h 



± 

L/R 



- ^L/R =•= *^L/fl- 



In the classical limit, 



neglecting the possibility of a finite bath polarization, 
^L/R^L/R ~ (^l/r)^ + e^L/i?)^- '^'^^ classical expres- 
sions for the hf-related terms in the effective Hamiltonian 
are then: 

- Uhl-h^R), 



Vss 

Vs^To 

Ha 

Vs9 



J 



Afi^^J^ 



H 



H 



h^) 



1, 1 



Mt 



Mt 



fJ-T 



J2 



)[(h 



±\2 



±n21 



iK) 



J 



4Mt- 

- — (M 

— ( — 



J2 



z • (h^ X hi) 



h-R)' 



veAh' [{hi 



±\2 






[{hi) 



(Cl) 
(C2) 

(C3) 

(C4) 

(hA)^] 
(C5) 

we can 
For ex- 



Replacing (/i|^ — /i|j)/2 by (Tz, and |h-'-| by a±, 
estimate the typical magnitude of these terms, 
ample, the S-Tq dephasing function due to Vss is given 
by 



Wisit) = 



<Phi 



(Phi 



'^^<^l.R 



exp [—i^t{h 



hi?] 



(C6) 



22 



Keeping in mind that a"^ — a'j^ j^ + a'^ j^, we calculate the 
Gaussian integrals and obtain the expression for PF55 (t) , 
which turns out to be equal to the one given in Eq. ([57l) . 
In other words, the semiclassical calculation of dephasing 
due to Vss is equivalent to neglecting the Larmor pre- 
cession of nuclear spins in the RDT calculation given in 
the main text. 

Similarly, the dephasing function due to Vse is given 
by 



where the index u labels the Wiger-Seitz (WS) unit cells, 
and fQ is the volume of the WS unit cell. The envelope 
wavefunction is normalized as /|^(r)p = i/q, and the 
function p{£^) parametrizes the "density of states" of hf 
couplings'^^ 



f J 



(D2) 



Wfeit) = 



with 



dh'r 



xWLih 



.lV 



dhl 






2TTaz,R 



As an example, if we assume that the envelope wavefunc- 
tion is a two-dimensional Gaussian: 



L,h^R)WRihi,h^^) 



(C7) 



WL{hl,h%) = 





■±-yLj 


y/iraL 




we have 





vi,(r) ^ ^^e-(-'+^')/2^'e(|z| - a/2) , (D3) 



xe 



-2tveAh'{(hlf + {hlf) 



(C8) 



The expression for Wr is analogous, only with the sign 
in front of ve inverted. After carrying out the integrals 
in Wl and Wr, we obtain 



WL/R{hl,h^j,) 



1 



l±4:iS,iAh'^aj ,„t 



(C9) 



'L/R- 



The decoherence function for a fixed value of Ah^ = 
26 can be obtained by multiplying the two expressions 
above. Since the coherence decay is faster when Ah^ is 
larger, we find the lower bound for the decay time by re- 
placing Ah^ with 2az- Assuming dots of equal sizes, we 
obtain the expression for Wgg{t) that is identical to the 
short-time (or single- isotope) result from Eq. (l66l) . 

In both calculations here we have neglected any non- 
trivial effects that the bath polarization can have on the 
free induction decay signal. These effects, discussed in 
Ref. [8l|, cannot be captured by the semiclassical approx- 
imation. However, at small nuclear polarizations consid- 
ered in this paper, they amount to a rather small quan- 
titative change in decay timescalCj^^ which justifies our 
approximation of neglecting them. 



^«)-tKi-«'- 



(D4) 



where 



m 



En l*(r«) 



27raL2 



(D5) 



As an example application, let us use the above p(^) to 
calculate the quantity appearing in Section IIII Al 



E !*(''") 



eV(Ode = ^ 



Nl, 



(D6) 



The prefactor 2 here depends on the exact shape of the 
wavefunction, but the ^N^^ scaling is a general result. 



Appendix E: Effective Dephasing Hamiltonians 
beyond the S-To Subspace 



Appendix D: The distribution of hyperfine couplings 
dependent on the wavefunction shape 

In the study of coupled electron-nuclear spin dynamics 
we often need to evaluate moments of the hyperfine in- 
teraction. Here we discuss the wave function dependence 
of these moments. Fr example, the nth moment can be 
written as a sum over the nuclear index i as 

i OL u 

= E"""^"/^>(Od?, (Dl) 



For single-spin qubits, exchange interaction is used to 
perform two-qubit gates and to transport spin states^i In 
these operations, generally all four two-spin states are in- 
volved. For simplicity, we also assume there is no inten- 
tional field gradient between the neighboring quantum 
dots (^T — 0). The presence of a field gradient would 
complicate gates such as a simple SWAP. Under these 
conditions, the effective dephasing Hamiltonians in the 
two-spin Hilbert space can be written as 



HsT = Hs\S){S\+Ht^\T+){T+\ 

+Ht„ \To) (ToI + Ht_ |r_) (r_| , (ei) 
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where 



Hs = -J + Ha- 



2(mt + J) 



Y,B^B,iriJ 



*j 



^W^,^"'^''''' 



-J + Ha + Vss + Hb 
-Mt4 

1 



^A-^E^^^.^r/; 



2(mt - J) 



T.B^B,iti- 



(E2) 



(E3) 



Ht^ - 






i?T 






2(/iT + J) 



Hb^b^i+i, 



2^^^^^^^^^^ 



i^A 



J 



(E4) 

(E5) 
(E6) 



This is the Hamiltonian we use to calculate the exchange 
gate errors due to hyperfine interaction with nuclear spins 
in Sec. Ivl 
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